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ABSTRACT 

We investigate the Ha and infrared star formation rate (SFR) diagnostics 
for galaxies in the Nearby Field Galaxy Survey (NFGS). For the 81 galaxies in 
our sample, we derive Ha fluxes (included here) from integrated spectra. There 
is a strong correlation between the ratio of far-infrared to optical luminosities 
L(FIR)/L(Ha) and the extinction E{B — V) measured with the Balmer decre- 
ment. Before reddening correction, the SFR(IR) and SFR(Ha) are related to 
each other by a power-law: SFR(IR) = (2.7 ± 0.3) SFR(Ha)^■^°^°■°^ Correc- 
tion of the SFR(Ha) for extinction using the Balmer decrement and a classi- 
cal reddening curve both reduces the scatter in the SFR(IR)-SFR(Ha) correla- 
tion and results in a much closer agreement between the two SFR indicators; 
SFR(IR) = (0.91 ± 0.04) SFR(Hacorr)^ °^^° °^. 

SFR(IR) and SFR(Ha) agree to ~ 10%. This SFR relationship spans 4 orders 
of magnitude and holds for all Hubble types with IRAS detections in the NFGS. A 
constant ratio between the SFR(IR) and SFR(Ha) for all Hubble types, including 
early types (SO-Sab), suggests that the IR emission in all of these objects results 
from a young stellar population. 

Subject headings: 
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1. Introduction 

Understanding the star formation history of the universe is the primary goal of much 
current research in astronomy (for example, Rowan-Robinson 2001; Cole et al. 2001; Baldry 
et al. 2002; Lanzetta et al. 2002; Rosa- Gonzalez, Terlevich, & Terlevich 2002). To obtain 
this understanding, a reliable estimate of the star formation rate in individual galaxies is 
required. Many calibrations of star-formation rate depend on the luminosity measured at 
various wavelengths including: radio (eg., Condon 1992; Cram et al. 1998; Haarsma et al. 
2000), IR (eg.. Hunter et al. 1986; Lchnert & Heckman 1996; Meurer et al. 1997; Kcnnicutt 
1998), optical (eg., Gallagher et al. 1989; Leitherer & Heckman 1995; Kennicutt 1998; Madau, 
Pozzetti, & Dickinson 1998; Rosa-Gonzalez, Terlevich, & Terlevich 2002), or UV (eg., Buat, 
Deharveng & Donas 1989; Deharveng et al. 1994; Leitherer, Robert, & Heckman 1995; Meurer 
et al. 1995; Cowie et al. 1997; Madau, Pozzetti, & Dickinson 1998; Rosa- Gonzalez, Terlevich, 
& Terlevich 2002). Unfortunately, the agreement among the SFR indicators at these different 
wavelengths is poor, and the underlying reasons for the differences are not well understood. 
To compute the SFR at any redshift and to relate SFRs at different redshifts, it is crucial to 
understand the physics underlying each SFR indicator and the source of the discrepancies 
among them. 

Here, we investigate two SFR indicators, the SFR(IR) and SFR(Hq;), for a large, ob- 
jectively selected sample of nearby galaxies. We aim to: (1) reduce the discrepancy between 
these two SFRs, and (2) obtain a better understanding of the relationship between the IR 
and Ha emission for all Hubble types. 

In star-forming galaxies. Ha photons are produced by gas ionized by young hot stars. 
The SFR(Hq;) is then a direct measure of the current SFR in galaxies, provided that red- 
dening is not significant. The SFR(Hq;) should agree with the SFR(IR) if the IR emission in 
star-forming galaxies arises from dust heated by the young hot stars. This scenario occurs in 
IR luminous starbursts (eg. Dopita et al. 2002), but the IR emission in normal star- forming 
galaxies appears more complex. 

The IR emission (l-lOOO/im) in normal star-forming galaxies may arise from three 
processes: (1) the emission from dust heated by young OB stars (Devereux & Young 1990; 
Devereux & Hameed 1997); (2) the emission from the photospheres or circumstellar envelopes 
of evolved stars such as red-giants undergoing mass-loss (Knapp, Gunn, & Wynn- Williams 
1992; Mazzei & dc Zotti 1994); (3) the "cirrus" emission from dust distributed throughout 
the optically thin, neutral interstellar medium, heated by the general stellar radiation field. 
The general stellar radiation field is composed of a combination of unabsorbed radiation 
from young OB stars inside the active star- forming regions (eg., Zurita, Rozas, h Beckman 
2000), and an underlying old stellar population consisting of A, F and G dwarfs, and K and 
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M giants (eg., Drapatz 1979; Helou 1986; Lonsdale & Helou 1987; Sauvage & Thuan 1994; 
Mazzei & de Zotti 1994). The relative contribution to the general stellar radiation field by 
the young and old stellar populations is unclear and may vary substantially from galaxy to 
galaxy depending on the details of the physical and dynamical environment and metallicity. 

The calibration of infrared luminosity as an indicator of the global star formation rate 
in a galaxy relies upon two assumptions; (a) that young stars dominate the radiation field 
throughout the UV to the visible (ie., that processes (2) and (3) above are negligible), and 
(b) that the dust opacity is infinite throughout a galaxy. If these assumptions hold, then 
the IR luminosity approximates the bolometric luminosity of the galaxy, and the SFR(IR) 
is a reliable estimate of the true SFR of the galaxy. 

For infrared-bright dusty star-forming galaxies, the young stellar population does pro- 
duce the dominant dust heating radiation field, and the dust opacity is high (Lonsdale & 
Helou 1987; Poggianti, Bressan, & Franceschini 2001; Dopita et al. 2002; Rosa- Gonzalez, 
Terlevich, & Terlevich 2002). Therefore, assumptions (a) and (b) above hold, and IR lumi- 
nosity is a more direct indicator of the true SFR than the Ha luminosity. Recently, Dopita et 
al. (2002) (hereafter D02) and Rosa- Gonzalez, Terlevich, & Terlevich (2002) (herafter R02) 
showed that using the B aimer decrement to correct the Ha luminosity for reddening results 
in better agreement between the SFR(IR) and SFR(Ha) . The D02 sample contains merging 
infrared starburst galaxies; the R02 sample consists of nearby young star-forming galaxies, 
most arc low metallicity late-type spirals or bhic compact dwarfs. It is unclear from these 
results whether the approximate agreement between SFR(IR) and SFR(Ha) after reddening 
correction also applies to a wider range of "normal" galaxies. 

The young (OB) star-forming population dominates the IR emission in late- type spirals. 
The nature of the IR emission in early-type spirals, on the other hand, remains controversial. 
Dust heating by the general stellar radiation field may make a significant contribution to 
the IR emission in early- types (eg., Lonsdale & Helou 1987; Sauvage & Thuan 1994; Mazzei 
& de Zotti 1994) bringing the use of SFR(IR) into question for these galaxies. Early-type 
galaxies may have small dust opacities further comphcating the use of SFR(IR). 

Kennicutt & Kent (1983) demonstrated that the Ha equivalent width increases system- 
atically from early-type to late-type galaxies. This increase suggests that early-type galaxies 
are deficient in young stars. Sauvage & Thuan (1994) provide additional evidence for this 
deficiency. They found that the B — V and U — B colors decrease and the ratio of H I 
mass to blue luminosity increases systematically from SO to Im type in the CfAl galaxy 
sample (Davis & Peebles 1983; Huchra et al. 1983). Furthermore, IR SED modeling of early- 
types appears to require a cooler component to the IR emission (Buat & Deharveng 1988; 
Rowan- Robinson & Crawford 1989), often attributed to the general stellar radiation field. 
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An intriguing theoretical argument by Inoue (2002) suggests that the SFR(IR) can be 
apphed to a wide range of galaxy types, including those with an older stellar population and 
small dust opacity. Inoue proposes a scenario in which the two effects of small dust opacity 
and a large cirrus contribution offset each other, making the SFR(IR) a good measure of the 
true SFR within a factor of 2 for a wide range of galaxy types. 

Although the general stellar radiation field may be important in many early-type spi- 
rals, H« images provide convincing evidence that at least some early-type galaxies contain 
significant recent star formation (Hameed & Devereux 1999). In addition, Tomita, Tomita & 
Saito (1996) and Dcvcrcux & Hamccd (1997) showed that the far infrarcd-to-blue luminosity 
ratio is independent of Hubble type for a large sample of nearby spirals. Their work supports 
the idea that early-type galaxies contain substantial recent star formation. 

A number of studies compare various SFR indicators for normal disk galaxies (eg., 
Kennicutt 1983; Hopkins et al. 2001; Chariot et al. 2002; Buat et al. 2002). Most of these 
authors cannot properly correct the Ha luminosity for reddening or stellar absorption because 
they lack integrated optical spectra with sufficient spatial coverage, signal-to-noise, and 
wavelength coverage. For example. Chariot et al. (2002) use had to assume an 'average' 
attenuation Ay = 1 for the galaxies in the Stromlo-APM survey where the slit-covering 
fraction is 40-50 percent. 

In this paper, we analyze integrated spectra from the Nearby Field Galaxy Survey 
(NFGS) (Jansen et al. 2000a), described in Section 2. We ensure that the integrated Ha 
luminosity is properly corrected for reddening and underlying stellar absorption. We compare 
uncorrected SFR(Hq;) with SFR(IR) in Section 3. Section 4 compares the corrected SFR(Ha) 
with the SFR(IR). The agreement between the two SFR indicators in the latter section is 
remarkable and is independent of Hubble type. We discuss the implications of this agreement 
in Section 5. 



2. Sample selection and cross-correlation with IRAS 

The NFGS is ideal for investigating relative IR and Ha star formation rates. The sample 
was selected objectively from the CfAl galaxy survey (Davis & Peebles 1983; Huchra et al. 
1983) with mB(o) < 14.5. The CfAl catalog is nearly complete within its selection limits 
and contains galaxies with a large range of absolute magnitude (—22 < < —13). Jansen 
et al. (2000a) chose a subsample of 198 nearby galaxies spanning the full range of Hubble 
type and absolute magnitude present in the CfAl catalog. To avoid a strict diameter limit, 
which might introduce a bias against the inclusion of low surface brightness galaxies in the 
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sample, they chose a radial velocity hmit, VLG(kms~"-'^) > IQ-^-^^-^ '^^^ (with respect to the 
Local Group standard of rest). To avoid a sampling bias favoring a cluster population, they 
excluded galaxies in the direction of the Virgo Cluster. The absolute magnitude distribution 
in the NFGS sample approximates the local galaxy luminosity function, while the distribution 
over Hubble type follows the changing mix of morphological types as a function of luminosity 
in the local galaxy population. Low surface brightness galaxies are overrepresented in the 
NFGS relative to the CfAl sample. 

Both integrated and nuclear spectrophotometry arc available for almost all galaxies in 
the NFGS sample, including reliably measured integrated Ha and H/? fluxes relative to the 
flux at 5500A (Jansen et al. 2000b; Jansen, Franx, & Fabricant 2001). The spectra have 
high enough S/N ratios to measure the H/3 fluxes even when E(B-V) is large. 

The integrated spectra typically cover 82± 7% of the galaxy. We calibrate the inte- 
grated fluxes, Fi.ei(A), to absolute fluxes, F(A), by careful comparison with B-band photom- 
etry. We converted total 5-magnitudes (extrapolated to infinite radius) to flux densities, 
F(4350A, phot), at the effective wavelength of the S-filter (4350 A). The uncertainties in 
the constant relating the magnitude to the flux of Vega, and in the effective wavelength of 
the S-filter for a given galaxy SED are ~2-3% combined. Because the position, size and 
orientation of the spectroscopic apertures are known, the i?-filtcr galaxy images allow us to 
compute the fractions, /, of the total i?-filtcr light that enters the spectroscopic apertures 
with 1.5-3% accuracy. The mean value of / = 0.82 ±0.07. The i?-filter flux density to which 
we must scale the spectra is then 

F(S,spec) = / X F(4350A,phot). 

Scaling the relatively flux calibrated spectrum, Frci(A), by the ratio of this flux density and 
that derived from the spectrum itself, we obtain the absolute flux calibrated spectrum: 

"'^^ - ■ (f„!;'a) t:(a)) 

where (Fj.ei(A) ® Tb(A)) denotes the mean relative flux density, derived by convolving the 
relative spectrum with the normalized B-flltcr response curve (Bessell 1979). The maximum 
error in the relative flux calibration is ~6.5%. The range in uncertainty in the total i?-flltcr 
magnitudes is ~1-10%, with a mean of 3.0% and rms of 2.7%. The errors in the absolute 
flux calibration of the spectrum are smallest at 4350A and increase toward both shorter and 
longer wavelengths. The expected error in the derived absolute Ha (H/3) emission hne fluxes 
is on average ~12% (~17%) and at most ~34% (~50%) for galaxies with EW(H/3) < -O.sA. 
The error in the H/3 flux is larger than that in Ha, because the H/3 emission hne is fainter 
and its measurement is complicated by the superposition of stellar absorpion. 
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As described in Jansen, Pranx, & Fabricant (2001), the line fluxes are corrected for 
Galactic extinction, using Burnstein & Hefles (1984; as listed in de Vaucouleurs et al. 1991). 
The H/3 emission is partially corrected for stellar absorption by ensuring that the limits 
of the flux measurement window are well inside the absorption trough, closely bracketing 
the emission line. We evaluated the residual absorption using spectra of galaxies without 
detectable emission. On average, we applied an additional correction of l.oA in equivalent 
width (EW) to H/3 and I.SA EW to Ha. 

Wc cross-correlated the NFGS with the IRAS Faint Source Catalogue (FSC) and Point 
Source Catalogue (FSC). The IRAS beam size at 60 /im is 1.5 arcminutes, and the IRAS 
positional uncertainty is 30 arcseconds. Because the NFGS positional uncertainty is only a 
few arcseconds and because all but a few of the NFGS galaxies are smaller than 3 arcmin 
in extent, we used a detection radius of 30 arcseconds. Any IRAS source associated with 
an NFGS object should be detected within this radius. We found 110 (56%) NFGS galaxies 
with IRAS sources within 30 arcseconds; 84 of these galaxies have moderate or good quality 
fluxes at 60 and 100//m, and Ho; and H/3 fluxes accurate to within 30%. 

To rule out the presence of AGN in the NFGS sample, we used the theoretical optical 
classiflcation scheme developed by Kewley et al. (2001). The optical diagnostic diagrams 
indicate that the global spectra of 81/84 NFGS galaxies are dominated by star formation. 
The spectra of the remaining 3 galaxies arc dominated by AGN: UGC 00545, Mrk 0205, and 
NGC 5940 host Seyfert I nuclei. These 81 galaxies constitute the sample we analyse here. 
We list the galaxies in Table 1. 



3. Infrcired versus Ha SFRs 

The development of SFR diagnostics has been an intense topic of research for over three 
decades (see Kennicutt 1998, for a review) (hereafter K98). Wc use the K98 SFR relations 
for IR and Ha because these relations arc used in many current SFR studies. We defer a 
more detailed analysis of other SFR diagnostics to a future paper. 

The IR SFR relation (K98) depends on the total infrared luminosity (8-1000//m) (Lir): 

SFR(IR)(MQyr-i) = 4.5 x 10"^^ Ljr (ergss-^). (1) 

Note that Lir, is referred to as LpiR in K98. For the remainder of this paper, we refer to 
the 40-120//m range as the far- infrared (FIR), and the l-1000//m range as the infrared (IR) 
regime. 
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The K98 SFR relation assumes that young stars dominate the radiation field throughout 
the UV-visible regimes, and that the optical depth is large. The infrared luminosity is then 
a good representation of the bolometric luminosity of the galaxy. Although this situation 
may hold for IR-luminous starburst galaxies, it does not necessarily apply to normal galaxies 
where the emission may be dominated by emission from a cooler, "infrared cirrus" component 
(eg., Lonsdale & Hclou 1987). 

In order to estimate the infrared star formation rate SFR(IR), we first calculate the 
far-infrared fiux (F(FIR): 40-120/im) from the IRAS 60 and 100/xm fiuxes (Fgo, Fioo) in units 
of Janskies (Helou et al. 1988): 

F(FIR) = 1.26 X 10-^^[2.58F6o + Fioo] (Wm'^) (2) 
The far-infrared luminosity L(FIR) is then calculated using: 

L(FIR) = 47rD2F(FIR) (3) 

where Dl is the cosmological "luminosity distance". We assume Ho=75 km s^^ Mpc^^ and 
qo = 0.5. For thermal energy distributions relevant for most galaxies, the FIR fiux and 
luminosity are essentially model independent. All galaxies in our sample have measured 60 
and lOO/xm fiuxes. 

New IR photometry and spectra from the Infrared Space Observatory (ISO) have lead 
to recent advances in infrared spectral energy distribution (SED) modelling (eg., Spinoglio, 
Paola & Malkan 2002; Dale ct al. 2001; Rowan-Robinson 2001; Calzctti ct al. 2000; Granato 
et al. 2000). It is now possible to obtain estimates of the total infrared luminosity from the 
IRAS 60 and 100//m fluxes (eg.. Chariot et al. 2002; Rosa-Gonzalez, Terlevich, & Terlevich 
2002). The uncertainties that arise from the conversion of 60 and 100//m IRAS fluxes to 
total (1-1000/xm) IR flux are likely to be modest in comparison with the uncertainties in the 
conversion from infrared luminosity to star formation rate (Gharlot et al. 2002). 

To transform L(FIR) to L(IR), we use the Galzetti et al. (2000) (hereafter GOO) for- 
mulation because it was developed for IR-luminous starburst galaxies and is therefore the 
most consistent relation to use in conjunction with the K98 SFR formula. GOO modelled the 
total-infrared flux for five low-redshift starburst galaxies with a two-component dust model. 
The dust model consists of two modified Planck functions of temperature T~ 40 — 55 K 
(warm dust) and T~ 20 — 23 K (cool dust), with an emissivity index e = 2. GOO finds that 
the ratio of the total dust IR emission (l-lOOO/^m) to the IRAS FIR emission is 



-8- 



L(IR) ~ 1.75 X L(FIR) 



(4) 



with little variation from galaxy to galaxy. Using this conversion, the K98 SFR(IR) becomes: 



Here, we define SFR(IR) as the SFR derived from the total (l-lOOO/xm) IR luminosity. 
We estimate the IR luminosity from the measured IRAS 60 and 100//m fluxes. We note that 
K98 defined the SFR(IR) over 8-1000//m. The COO conversion apphes to the range l-1000//m. 
The contribution to the total infrared luminosity from the 1-8/xm regime is expected to be 
of the order of a few percent (Calzetti et al. 2000; Dale et al. 2001). 

The Ha SFR relation (K98) comes from evolutionary synthesis models, assuming solar 
metallicity and no dust. The total integrated stellar luminosity shortward of the Lyman 
limit is re-emitted in the nebular emission lines. The Ha SFR relation is: 



For most star-forming galaxies, dust absorption leads to the underestimation of the SFR 
derived from the Ha luminosity (eg., Lonsdale & Helou 1987). 

Figure 1 shows the SFRs derived using the K98 IR and Ha SFR relations for the NFGS 
galaxies. There is obviously a strong correlation between the two SFRs. Upper hmits for 
those objects with measurable Ha fiuxes but without IRAS detections at 60 and 100//m are 
consistent with this correlation. 

We fit a straight fine to the logarithm of the SFRs for the objects with IRAS detections 
using the numerical recipes fitexy routine in IDL. This routine uses linear least-squares 
minimization and includes error estimates for both variables. We assumed errors of ~60% 
for the non-reddening corrected SFR(Ha) and ~30% for SFR(IR). The large ~60% error 
allows for the systematic error introduced by failure to correct for reddening. The resulting 
fit (dotted line in Figure 1) has the form: 



SFR(IR) (Moyr-i) 



4.5 X 10-4^L(1R) 

7.9 X 10~^^L(FIR) (ergss^^) 



(5) 



SFR(Ha)(MQyr-^) = 7.9 x 10-^^L{lla) (ergss"^) 



(6) 



log[SFR(IR)] = (1.30 ± 0.06) log[SFR(Ha)] + 0.43 ± 0.05 



(7) 



or: 



SFR(IR) = (2.7 ± 0.3)SFR(Ha) 



1.30±0.06 



(8) 



-9- 



The rms distance perpendicular to this line (in the log) gives a measure of the degree of 
scatter. The rms dispersion in Figure 1 is 0.17 in the log. The Spearman Rank correlation test 
for SFR(IR) and SFR(Hq;) gives a correlation coefficient of 0.95. The two-sided probability 
of finding a value of 0.95 by chance is formally 0.0 (< 10~^°), confirming the very strong 
correlation between the two SFR indicators. 

We plot a histogram of the logarithm of the ratio of the two SFRs, log[SFR(IR)/SFR(Ha)], 
in Figure 2. Markers at the top of Figure 2 show the means of the early and late type SFR 
ratio distributions, 0.50 ± 0.07 and 0.34 ± 0.04, respectively. We note that the mean for the 
early types is larger than the mean for the late types. 

4. SFRs and Reddening 

Figures 1 and 2, and equation 8 demonstrate that the SFR(IR) estimate exceeds SFR(Ha) 
by a factor of ~ 3, as observed by Chariot et al. (2002). The ratio between SFR(IR) and 
SFR(Ha) increases with SFR. If the SFR(IR) and SFR(Ha) are both valid measures of the 
true SFRs, then a physical process not yet taken into account must bias one relative to 
the other. This process must also cause the bias to increase at higher SFRs. One obvious 
candidate is reddening. We test this conjecture by measuring the difference in the extinction 
at two wavebands, E{B — V). 

We calculated E{B — V) in the emission-line gas for the 81 IRAS-detected galaxies 
in the NFGS sample for which flux measurements for both Ha and H/3 are available. We 
assume an intrinsic Hq;/H/3 ratio of 2.85 (case B recombination at T= lO^K and electron 
density ~ 10^ — lO^cm"^; Osterbrock (1989)) and used the Whitford reddening curve as 
parameterized by Miller & Mathews (1972). 

Figure 3 shows the relationship between the ratio of the two luminosities L(FIR)/L(Hq;) 
and E{B — V). The Spearman Rank correlation correlation coefficient is 0.86. The two-sided 
probability of obtaining a value of 0.86 by chance is ~ 1.2 x 10~^^, supporting the very strong 
correlation between L(FIR)/L(Hq;) and E{B - V). 

The fitexy least squares line fitting routine gives 

= (0.62 ± 0.08) log[E(B - V)] + 2.66 ± 0.06. (9) 

This relationship is model independent. If we use E{B — V) to correct L(Hq;) for extinction, 
the difference between SFR(IR) and SFR(Hq;) could be reduced, supporting similar results by 



log 



L(FIR) 
L{Ra) 
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Dopita et al. (2002) and Rosa- Gonzalez, Terlevich, & Terlevich (2002) for young starforming 
galaxies. Infact, after making the correction, we do arrive at much better agreement between 
the corrected Ha star formation rate (SFR(HQ;corr)) and the IR star formation rate (Figure 4). 
This relation is striking because many nearby field galaxies are cool IRAS sources. Only 6 
out of the 81 galaxies in Figure 4 have a 60;um to lOOyum flux ratio greater than 0.6, typical 
for objects dominated by hot stars (Bothun, Lonsdale, & Rice 1989). 

We fit a straight line to the NFGS data in Figure 4 using the same method as for 
Figure 1. The errors are ~ 30% for both SFR(Hq;) and SFR(IR). The resulting fit is 

log[SFR(IR)] = (1.07 ± 0.03) log[SFR(HQ;,orr)] - 0.04 ± 0.02 (10) 

or, 

SFR(IR) - (0.91 ± 0.04)SFR(Ha;corr)^-°^^° °^ (11) 

The Spearman rank correlation coefficient for (SFR(HQ;corr)) and SFR(IR) is 0.98, with 
the 2-sided probability of obtaining this value of < 10"^°. The rms scatter of the data 
around this line decreases to 0.12 in the log. Clearly correcting Ha for reddening results in 
a closer correspondence between SFR(IR) and SFR(Ha), and reduces the scatter observed 
in Figure 1. 

If the general stellar radiation field contributes significantly to the IR emission in early- 
types, we would expect the late-types to display a larger scatter than early types. We find 
the contrary: the rms scatter for the early-types (SO-Sab) is 0.08, compared to 0.13 for the 
late- types (Sb-Im). 

We plot a histogram of the logarithm of the ratio of the two SFRs, log(SFR(IR) /SFR(Ha)) , 
in Figure 5. After correction for reddening, the means of the two SFR ratios for early and 
late types are approximately equal: —0.02 ± 0.03 and —0.07 ± 0.02 respectively. Greater 
reddening for the early types appears to explain the offset in Figure 2. The overall mean in 
the ratios is -0.05 ± 0.02. 

By selecting early-types with infrared detections, we bias our early-type sample towards 
galaxies with more recent star-formation than is normal in these galaxies. The number of 
early-types with 60 and 100 fim IRAS detections is 17/59 (29%), compared with the number 
of late-types with detections 73/112 (65%). This detection rate supports previous studies 
which showed that early- types in general have less recent star formation than late types (eg., 
Kennicutt & Kent 1983; Sauvage & Thuan 1994; Bendo et al. 2002). The strong agreement 
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and very small scatter between SFR(IR) and SFR(Hq;) for the infrared early-types implies 
that they are dominated in both Ha and the IR by recent star formation. 

We obtain SFR estimates by multiplying the luminosity at the relevant wavelength by a 
constant. Thus, the SFR correlation (Figure 4) is really equivalent to a correlation between 
the far-infrared and Ha luminosities. The line of best fit to the luminosity correlation is; 



log 



L(FIR) 



^0 



;i.07±0.02)log 



L(Hae 



^0 



+ (1.4 dz 0.2) 



(12) 



or. 



L(FIR) 



(25± 12) 



-■0 



L(Ha, 



^0 



(1.07±0.02) 



(13) 



For the K98 star formation rate indicators to produce equivalent IR and Ha star for- 
mation rates, the ratio L(FIR)/L(Hacorr) multiplied by the ratio of the constants in the K98 
relations (including the COO factor of 1.75 to convert L(FIR) to L(IR)) should be ~ 1. The 
mean L(FIR)/L(Hacorr) for our sample is 96 ± 4. Multiplying the mean L(FIR)/L(Hacorr) 
by ratio of the K98 relations gives (96 ±4) (7.9 x 10-^^)/(7.9 x lO'^^) ^ 0.96 ±0.04. Clearly, 
the combination of the COO conversion and the ratio of the K98 coefficients gives an excellent 
estimate of the relationship between L(FIR)/L(Hacorr) and SFR(IR)/SFR(Hacorr)- 

It is interesting that such a strong FIR-Ha correlation occurs. In the simplest pictTu^c. 
the correlation indicates that a single physical process is responsible for both the FIR and 
Ha emission. In reality, the FIR and Ha radiation may be produced in quite different 
environments. The Ha radiation comes from gas ionized by the UV radiation from young, 
hot OB stars. This process occurs within many H II regions throughout a galaxy. The FIR 
radiation source is more complex. 

The far-IR spectra of galaxies result from large dust grains of diameter 0.01-0.3. The 
spectrum from these grains is characteristic of thermal emission at 15-20K (sec Draine 2001, 
for a review), and accounts for nearly all the emission at A > 60/im. Some authors believe 
that an older stellar population can make a significant contribution to the far-IR emission 
(eg., Lonsdale & Helou 1987; Sauvage & Thuan 1994; Mazzei & de Zotti 1994) for early-type 
galaxies. However, our Ha-FIR luminosity correlation implies that the dominant UV and 
optical heating source for this emission is the young stellar population. This hypothesis is in 
agreement with various other physical correlations linking the FIR emission to young stellar 
populations, such as the radio-FIR correlation (eg., Gavazzi, Cocito, & Vettolani 1986), and 
the 12yum-FIR luminosity correlation (Shapley, Fabbiano, & Eskridge 2001). Both of these 
correlations hold for all galaxy types. 
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These correlations suggest that the gas and dust must be closely coupled in all FIR- 

bright galaxies. Malhotra et al. (2001) found a strong correlation between L[oi]/L[cii] and 
F(60/im)/F(100/im), indicating that both gas and dust temperatures increase together. The 
most plausible explanation for this effect is that the majority of the dust and gas heating in 
galaxies occurs very close to the actively star-forming regions. 

Modeling by Popescu &: Tuffs (2002) for normal star-forming galaxies strengthens this 
explanation. The Popescu et al. models predict that radiation from H II regions is the dom- 
inant energy source for the dust emitting at 60/im. Even at lOOyum, H II regions contribute 
~ 1/3 of the radiation source. The other 2/3 is produced by the diffuse UV and optical 
radiation, which could also contain some contribution from the young stellar population. 

Even if the FIR emission is dominated by the young stellar population, one would 
expect the ratio of gas and dust masses to vary from one H II region to another, and from 
one galaxy to another. Similarly, the dust geometry and composition must vary. If these 
variations occur, they do not appear to have an effect on the relationship between the global 
FIR and Ho; luminosities. 



5. Conclusions 

We analyze a sample of 81 far-infrared detected galaxies from the Nearby Field Galaxy 
Survey. The NFGS is an objectively selected sample of galaxies in the CfAl catalog. The 
sample contains a representative mix of all galaxy types and morphologies represented in the 
CfAl survey. In addition, the NFGS is unique in the quahty and coverage of the integrated 
spectra. We have taken extensive care to ensure that the integrated Ho; flux for each galaxy is 
as reliable as possible. We have corrected these fluxes carefully for extinction and underlying 
stellar absorption. 

To compare the IR and Ha star formation rates, we make the following standard as- 
sumptions: (1) the Whitford reddening law is appropriate for all galaxies in our sample, (2) 
the star formation rates are proportional to the relevant luminosities, and (3) that the IR 
SED of the galaxies in our sample can be approximated by the Calzetti et al. (2000) relation. 

With the assumed reddening, the data provide a direct measurement of the ratio of 
the IR and Ha SFR constants which transform the luminosities to star formation rates. 
Correction for reddening brings the IR and Ha SFRs into agreement to within 10% on 
average, and reduces the scatter. 

Correction for reddening also brings the mean of SFR(IR)/SFR(Ha) for early and late 
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types into agreement. One average, early types in this sample are more heavily reddened 
than late types. It is intriguing that the ratio between the SFR(IR) and SFR(Hq;) is the same 
for early and late types. The IR emission from early-type spirals may contain a significant 
contribution from the general stellar radiation field. This contribution should increase the 
observed IR luminosity. Inoue (2002) suggested that this effect is offset by a lower dust 
opacity. This compensation would require a remarkable conspiracy to hold for all galaxy 
types. The physical properties of dust and the nature of the general stellar radiation field 
probably vary from galaxy to galaxy. The variation in the general stellar radiation field 
should result in a larger dispersion around the SFR relation for early types than late-types. 
We do not observe this variation. Rather, the infrared NFGS early-type galaxies have a 
smaller dispersion than the late-types. 

Previous studies have shown that early-type galaxies contain a lower fraction of young- 
to-old stellar populations than late-types (eg., Kennicutt & Kent 1983; Sauvage & Thuan 
1994; Davis & Peebles 1983; Huchra et al. 1983). This result docs not determine what combi- 
nation of young-to-old populations produces the IR emission in each galaxy. The agreement 
of SFR(IR) with SFR(Hq;) (or L(FIR) with L(Hq;)) for both early types and late-types im- 
plies that the FIR emission results from the same young stellar population which produces 
the Ha emission. This hypothesis is supported by other physical correlations linking the 
FIR emission to the young stellar population, including the radio-FIR correlation (Gavazzi, 
Cocito, & Vettolani 1986) and the 12//m-FIR luminosity correlation (Shapley, Fabbiano, & 
Eskridge 2001). These correlations suggest that the dust is heated in close proximity to 
the young star-forming regions. In addition, the FIR-Ha luminosity correlation implies that 
the relationship between the global dust and gas properties is universal for all FIR-detected 
galaxies in the NFGS. 

Further observations and modeling are required to prove that the young stellar popu- 
lation is responsible for FIR emission in all galaxy types, rather than a compensation effect 
as suggested by Inoue (2002). Near-infrared observations would enable separation of the 
young and old stellar populations. The Space Infrared Telescope Facility (SIRTF) will help 
determine whether the global dust and gas properties are truly universal for all FIR-detected 
galaxies in the NFGS. 

We enjoyed discussions about many aspects of this paper with Rob Kennicutt, Dan 
Fabricant, and Stcphanc Chariot. L. J. Kcwlcy is supported by a Harvard-Smithsonian CfA 
Fellowship. M. J. Geller is supported by the Smithsonian Institution. 
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Ha uncorrected for reddening 




0.01 0.10 1.00 10.00 100.00 

SFR (Ha) 

Fig. 1. — A comparison of the SFRs derived using the infrared and Ha relations in Kennicutt 
(1998). The strong correlation spans 4 orders of magnitude. The solid line is y=x, and shows 
where the data would lie if both SFR indicators agreed. The Ha SFR has not been corrected 
for extinction. The dashed line is the best-fit to all of the data. The legend indicates the 
Hubble type. 
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-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 
log[ SFR(IR) / SFR(Ha) ] uncorrected 



1.0 1.2 



Fig. 2. — The ratio of the infrared to Ha star formation rates. The Ha flux has not been 
corrected for reddening. The two vertical lines at the top of the Figure indicate the mean of 
the early and late-type distributions. 
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Fig. 3. — The relationship between the reddening E{B — V) derived from the Balmer decre- 
ment F(Hci;)/F(H/3) and the ratio of the IR to Ha luminosities. 
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Ha corrected for reddening 




0.01 0.10 1.00 10.00 100.00 

SFR (Ha corrected) 

Fig. 4. — As in Figure 1, with the Ha flux corrected for reddening using the Whitford 
reddening curve as parameterized by Miller & Mathews (1972). The SFR(IR) and SFR(Hq;) 
are nearly equal. The legend indicates the Hubble types. The galaxies with the largest SFRs 
are late-type spirals. 
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■1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 
log[ SFR(IR) / SFR(Ha) ] corrected 



0.8 1.0 



Fig. 5. — The ratio of the infrared to Ha star formation rates. The Ha flux has been 
corrected for reddening. The two vertical lines at the top of the Figure indicate the mean of 
the early and late-type distributions. 
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Table 1. The NFGS Sample with 60 and 100/im IRAS detections 



ID 


Name 


cz 


Type 


E(B-V) 




r 

Floo'' log ^ 


Lcorr(Ha)l , 


rL(FIR)l 


SFR 


SFR 


SFR 


Lq J 








km/s 






Jy 


Jy 


Ha 


Hctcorr 


IR 


4 


UGC 00439 


5302 


Sa 


0.44 


1.14 


2.80 


8.13 


10.04 


1.50 


4.13 


3.40 


5 


UGC 00484 


4859 


Sb 


0.37 


0.91 


2.50 


7.94 


9.90 


1.15 


2.69 


2.42 


16 


UGC 01154 


7756 


Sbc 


0.43 


0.56 


1.96 


8.27 


10.15 


2.14 


5.72 


4.30 


17 


UGC 01155 


3158 


Sbc 


0.30 


0.49 


1.21 


7.35 


9.23 


0.35 


0.69 


0.52 


19 


NGC 695 


9705 


Sc 


0.75 


7.69 


12.84 


9.41 


11.33 


14.27 


79.67 


65.02 


21 


UGC 01551 


2669 


Sdm 


0.28 


0.48 


1.50 


7.42 


9.13 


0.42 


0.81 


0.41 


23 


IC 197 


6332 


Sbc 


0.42 


0.97 


2.31 


8.09 


10.12 


1.43 


3.79 


4.08 


25 


UGC 01630 


4405 


Sb 


0.44 


1.55 


3.26 


7.91 


9.99 


0.91 


2.50 


2.98 


27 


NGC 927 


8258 


Sc 


0.54 


0.82 


2.29 


8.52 


10.31 


2.98 


10.26 


6.32 


28 


UGC 01945 


1762 


Sdm 


0.24 


0.24 


0.79 


6.53 


8.48 


0.06 


0.10 


0.09 


37 


UGC 04713 


9036 


Sb 


0.96 


0.61 


2.81 


8.80 


10.39 


2.11 


19.37 


7.56 


39 


NGC 2780 


1951 


Sab 


0.34 


0.36 


1.04 


6.63 


8.71 


0.06 


0.13 


0.16 


43 


NGC 2844 


1486 


Sa 


0.37 


0.48 


1.46 


6.66 


8.61 


0.06 


0.14 


0.13 


44 


NGC 3011 


1517 


so/a 


0.09 


0.23 


0.61 


6.23 


8.28 


0.04 


0.05 


0.06 


45 


NGC 3009 


4666 


Sc 


0.30 


0.27 


0.82 


7.36 


9.35 


0.35 


0.71 


0.69 


46 


IC 2520 


1226 


../Pec 


0.46 


3.48 


6.73 


7.19 


9.21 


0.17 


0.48 


0.50 


47 


UGC 05354 


1172 


Sm 


0.04 


0.43 


0.92 


6.57 


8.29 


0.11 


0.11 


0.06 


48 


NGC 3075 


3566 


Sc 


0.15 


0.79 


1.72 


7.36 


9.52 


0.50 


0.71 


1.01 


49 


UGC 05378 


4185 


Sb 


0.46 


0.48 


1.30 


7.60 


9.48 


0.43 


1.24 


0.94 


50 


NGC 3104 


604 


Im 


0.03 


0.37 


1.00 


6.12 


7.69 


0.04 


0.04 


0.02 


53 


UGC 05522 


1228 


Sc 


0.08 


0.26 


0.85 


6.44 


8.20 


0.07 


0.08 


0.05 


56 


NGC 3213 


1412 


Sbc 


0.40 


0.25 


0.79 


6.52 


8.29 


0.04 


0.10 


0.06 


57 


NGC 3264 


929 


Sdm 


0.00 


0.46 


1.02 


6.49 


8.12 


0.11 


0.11 


0.04 


58 


NGC 3279 


1422 


Sc 


0.54 


1.83 


5.37 


6.90 


9.15 


0.07 


0.25 


0.43 


59 


UGC 05744 


3338 


Sc 


0.29 


0.86 


1.48 


7.30 


9.44 


0.32 


0.62 


0.84 


60 


UGC 05760 


2997 


Scd 


0.43 


1.05 


3.09 


7.42 


9.55 


0.30 


0.82 


1.10 


61 


IC 2591 


6755 


Sbc 


0.20 


0.59 


1.58 


7.94 


9.99 


1.68 


2.67 


2.99 


63 


UGC 05798 


1534 


Sc 


0.13 


0.36 


0.62 


6.42 


8.40 


0.06 


0.08 


0.08 


65 


CGCG 212-039 


10693 


Sbc 


0.64 


0.43 


1.32 


7.88 


10.28 


0.53 


2.33 


5.88 


73 


NGC 3510 


704 


Sd 


0.05 


0.60 


1.46 


6.38 


8.01 


0.07 


0.07 


0.03 


76 


IRAS F11041+5127 


2204 


Sc 


0.06 


0.44 


1.49 


6.91 


8.94 


0.22 


0.25 


0.27 


77 


IC 673 


3851 


Sa 


0.37 


0.59 


1.90 


7.57 


9.54 


0.49 


1.15 


1.07 


82 


NGC 3633 


2553 


Sa 


0.82 


3.08 


5.34 


7.75 


9.78 


0.26 


1.72 


1.83 


88 


UGC 06575 


1225 


Sc 


0.11 


0.29 


0.76 


6.16 


8.19 


0.03 


0.04 


0.05 


89 


NGC 3795 


1154 


Scd 


0.00 


0.23 


0.68 


6.08 


8.07 


0.05 


0.05 


0.04 


90 


UGC 06625 


10964 


Sc 


0.28 


1.23 


3.39 


8.56 


10.73 


5.82 


11.06 


16.59 


91 


NGC 3795A 


1091 


Sc 


0.27 


0.40 


1.34 


6.21 


8.29 


0.03 


0.05 


0.06 


96 


UGC 06805 


1033 


SO 


0.23 


0.52 


0.73 


6.29 


8.18 


0.04 


0.06 


0.05 


98 


IC 746 


5027 


Sb 


0.30 


0.50 


1.13 


7.61 


9.62 


0.63 


1.25 


1.29 


100 


NGC 3978 


9978 


Sbc 


0.38 


2.24 


6.88 


8.84 


10.94 


8.77 


21.27 


26.61 


105 


UGC 07020A 


1447 


SO 


0.31 


1.82 


2.29 


7.06 


9.00 


0.17 


0.35 


0.31 


107 


NGC 4120 


2251 


Sc 


0.14 


0.47 


1.58 


6.98 


8.99 


0.21 


0.30 


0.30 


109 


NGC 4141 


1980 


Sc 


0.07 


0.31 


0.85 


7.22 


8.64 


0.43 


0.51 


0.14 


110 


NGC 4159 


1761 


Sdm 


0.29 


0.81 


1.56 


7.02 


8.89 


0.17 


0.32 


0.24 


112 


NGC 4238 


2771 


Sc 


0.14 


0.49 


1.08 


7.23 


9.09 


0.38 


0.52 


0.38 


114 


UGC 07358 


3639 


Sc 


0.32 


0.46 


1.55 


7.40 


9.40 


0.37 


0.78 


0.77 


122 


UGC 07690 


540 


Sdm 


0.01 


0.54 


1.21 


5.88 


7.72 


0.02 


0.02 


0.02 


124 


NGC 4509 


907 


Sm 


0.04 


0.61 


0.52 


6.45 


8.07 


0.08 


0.09 


0.04 
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Table 1 — Continued 



ID 


Name 


km/s 


Tvne 


E(B-V) 


^ bU 

Jy 


^ iUU -'■'-'& 

Jy 


Lcorr(Hci)l 


. rL(FIR)l 


SFR 
Ha 


SFR 

Hctcorr 


SFR 
IR 


L0 J 


L0 J 


125 


UGC 07761 


6959 


Sb 


0.43 


0.96 


2.50 


8.14 


10.22 


1.57 


4.24 


5.09 


127 


NGC 4758 


1244 


Sbc 


0.39 


0.95 


2.53 


6.79 


8.72 


0.08 


0.19 


0.16 


133 


UGC 08231 


2460 


Sb 


0.07 


0.36 


0.97 


6.96 


8.90 


0.24 


0.28 


0.24 


134 


IC 4213 


815 


Scd 


0.12 


0.29 


0.71 


6.08 


7.82 


0.03 


0.04 


0.02 


135 


UGC 


08400 


3396 


Scd 


0.27 


0.42 


1.22 


7.26 


9.26 


0.30 


0.56 


0.56 


140 


NGC 


5230 


6855 


Sc 


0.33 


1.56 


4.71 


8.45 


10.45 


4.06 


8.61 


8.67 


141 


UGC 


08630 


2364 


Sm 


0.16 


0.60 


1.44 


6.98 


9.06 


0.20 


0.29 


0.35 


142 


NGC 


5267 


5941 


Sb 


0.91 


0.22 


1.12 


7.89 


9.61 


0.30 


2.41 


1.25 


144 


NGC 


5338 


777 


SO 


0.75 


0.39 


0.63 


6.09 


7.83 


0.01 


0.04 


0.02 


145 


NGC 


5356 


1397 


Sb 


0.99 


0.54 


2.40 


7.35 


8.71 


0.07 


0.69 


0.16 


148 


NGC 


5425 


2062 


Sc 


0.24 


0.37 


1.40 


6.89 


8.84 


0.14 


0.24 


0.21 


151 


NGC 


5491 


5845 


Sc 


0.46 


0.68 


2.19 


7.98 


9.97 


1.00 


2.91 


2.84 


153 


NGC 


5541 


7698 


Sc 


0.44 


2.63 


5.12 


8.52 


10.69 


3.71 


10.30 


14.91 


156 


UGC 


09356 


2234 


Sc 


0.12 


0.40 


0.97 


6.89 


8.84 


0.18 


0.24 


0.21 


158 


NGC 


5762 


1788 


Sa 


0.10 


0.20 


0.74 


6.59 


8.44 


0.10 


0.12 


0.08 


159 


UGC 


09560 


1215 


Pec 


0.00 


0.71 


1.24 


6.64 


8.49 


0.15 


0.15 


0.10 


162 


UGC 


09660 


626 


Sc 


0.16 


0.26 


0.61 


5.64 


7.54 


0.01 


0.01 


0.01 


164 


IC 1100 


6561 


Scd 


0.47 


1.03 


3.18 


8.07 


10.24 


1.23 


3.63 


5.30 


165 


NGC 5874 


3128 


Sbc 


0.48 


0.47 


1.85 


7.46 


9.31 


0.29 


0.88 


0.63 


166 


NGC 5875A 


2470 


Sc 


0.12 


0.51 


1.11 


6.87 


9.01 


0.17 


0.23 


0.31 


168 


IC 1124 


5242 


Sab 


0.60 


0.95 


2.71 


7.92 


9.99 


0.65 


2.58 


2.99 


170 


UGC 09896 


6461 


Sc 


0.21 


0.28 


1.17 


7.71 


9.74 


0.98 


1.60 


1.68 


172 


IC 1141 


4458 


so/a 


0.58 


1.34 


2.00 


7.65 


9.87 


0.36 


1.36 


2.29 


174 


NGC 


6007 


10548 


Sbc 


0.31 


0.70 


2.03 


8.40 


10.47 


3.84 


7.78 


8.98 


175 


UGC 


10086 


2191 


Sc 


0.25 


0.64 


1.28 


6.91 


8.99 


0.14 


0.25 


0.30 


179 


NGC 


6131 


5054 


Sc 


0.12 


0.72 


2.42 


7.51 


9.87 


0.74 


0.99 


2.30 


184 


NGC 


7328 


2827 


Sab 


0.64 


1.16 


3.95 


7.65 


9.58 


0.31 


1.38 


1.17 


186 


UGC 


12178 


1925 


Sdm 


0.15 


1.08 


2.72 


7.11 


9.15 


0.28 


0.40 


0.43 


187 


UGC 


12265 


5682 


SO 


0.57 


1.62 


2.53 


8.02 


10.17 


0.87 


3.23 


4.57 


189 


NGC 


7460 


3296 


Sb 


0.44 


1.82 


4.98 


7.93 


9.86 


0.94 


2.61 


2.22 


192 


UGC 


12519 


4380 


Sd 


0.45 


0.76 


2.59 


7.86 


9.78 


0.79 


2.25 


1.85 


193 


NGC 


7620 


9565 


Scd 


0.52 


2.41 


5.74 


9.01 


10.88 


9.35 


31.29 


23.13 


195 


IC 1504 


6306 


Sb 


0.85 


1.40 


4.12 


8.56 


10.32 


1.58 


11.23 


6.48 



*IRAS Fluxes at 60 and 100/im in units of Janskies 



